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£H \ Abstract 

• We consider the class of protocols that can be implemented by local quantum operations 

| and classical communication (LOCC) between two parties. In particular, we focus on the task 

. of discriminating a known set of quantum states by LOCC. Building on the work in the paper 

Quantum nonlocality without entanglement [BDF + 99], we provide a framework for bounding 
1 ^ 1 . the amount of nonlocality in a given set of bipartite quantum states in terms of a lower bound 

on the probability of error in any LOCC discrimination protocol. We apply our framework 
to an orthonormal product basis known as the domino states and obtain an alternative and 
simplified proof that quantifies its nonlocality. We generalize this result for similar bases in larger 
dimensions, as well as the "rotated" domino states, resolving a long-standing open question 
^ [BDF+99]. 
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1 Introduction 



The 1999 paper Quantum nonlocality without entanglement [BDF+99] exhibits an orthonormal basis 
S C C 3 (g>C 3 of product states, known as domino states, shared between two separated parties. When 
the parties are restricted to perform only local quantum operations and classical communication 
(LOCC), it is impossible to discriminate the domino states arbitrarily well [BDF + 99]. In such cases 
we say that perfect discrimination cannot be achieved with asymptotic LOCC. Moreover, [BDF + 99] 
also quantifies the extent to which any LOCC protocol falls short of perfect discrimination of the 
domino states. 

This result spurred interest in state discrimination with LOCC. Several alternative proofs 
[WH02, GV01, Coh07] of the impossibility of perfect LOCC discrimination of the domino states 
were given along with many other results concerning perfect state discrimination (e.g., [BDM + 99, 
WSHVOO, GKR+01, GV01, VSPM01, CY01, CY02, WH02, DMS+03, CL03, HSSH03, HM03, 
Fan04, GKRS04, Che04, CL04, JCY05, Wat05, Nat05, NC06, DFJY07, FS09, DFXY09, DXY10]). 
However, the problem of asymptotic LOCC state discrimination has not received much attention 
since the initial study of nonlocality without entanglement [BDF + 99]. 

The main motivation for our work is to better understand the phenomenon of quantum nonlo- 
cality without entanglement. More concretely, our goals are to 

• simplify the original proof, 

• render the technique applicable to a wider class of sets of bipartite states, 

• exhibit new classes of product bases that cannot be asymptotically (as opposed to just per- 
fectly) discriminated with LOCC, 

• pin down where exactly the difference between LOCC and separable operations lies, and 

• investigate the possibility of larger gaps between the sets of LOCC and separable operations. 

In particular, we seek to exhibit quantitative gaps between the classes of LOCC and separable 
operations. Separable operations often serve as a relaxation of LOCC operations and such gaps 
show how imprecise this relaxation can be. The rationale behind this relaxation is that separable 
operations have a clean mathematical description whereas LOCC operations can be much harder 
to understand. 

There is also an operational motivation to quantify the difference between separable measure- 
ments and those implemented by asymptotic LOCC: the former are precisely the measurements 
that cannot generate entangled states, while the latter are those that do not require entangle- 
ment to implement [BDF + 99, KTYI07, Koa09]. Thus, a separable measurement that cannot be 
implemented by asymptotic LOCC uses entanglement irreversibly. 

Our contributions 

In this paper, we develop a framework for obtaining quantitative results on the hardness of quantum 
state discrimination by LOCC. More precisely, we provide a method for proving a lower bound on 
the error probability of any LOCC measurement for discriminating states from a given set S. 

Our first main contribution (Theorem 2) is that any LOCC measurement for discriminating 
states from a set S errs with probability Perror 

— J?|S1 F ' wri ere r\ is a constant that depends on S 
(see Definition 3.4). Intuitively, r\ measures the nonlocality of S. 
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Our second main contribution is a systematic method for bounding the nonlocality constant 
T) for a large class of product bases. Together with the above theorem, this lets us quantify the 
hardness of LOCC discrimination for the following bases of product states: 

1. domino states, the original set of nine states in 3 x 3 dimensions first considered in [BDF + 99], 

have terror 

> 1.9 x 10~ 8 ; 

2. domino-type states, a generalization of domino states to higher dimensions corresponding to 
tilings of a rectangular dA x oIb grid by tiles of size at most two, have terror 

where D is a property of the tiling that we call "diameter" ; 

3. 0-rotated domino states, a 1-parameter family that includes the domino states and the stan- 
dard basis as extreme cases, have p error > 2.4 x 10~ n sin 2 29 (determining whether these 
states can be discriminated perfectly by LOCC and finding a lower bound on the probability 
of error were left as open problems in [BDF + 99]). 

The rest of the paper is organized as follows. In Section 2 we introduce notation, give back- 
ground on LOCC measurements and state discrimination, and summarize related prior work. In 
Section 3 we introduce our general framework for lower bounding the error probability of LOCC 
measurements, and in Section 3.5 we prove Theorem 2. In Section 4 we consider the case where S is 
a product basis and propose a method for bounding the nonlocality constant rj by another quantity 
that we call "rigidity." Our approach is based on a description of sets of bipartite states in terms 
of tilings. In Section 5 we define the three classes of states mentioned above and prove a bound on 
the rigidity of the domino states; bounds on the rigidity of the domino-type states and the rotated 
domino states appear in Appendices A and B, respectively. Finally, we discuss limitations of our 
framework in Section 6 and conclude with a discussion of open problems in Section 7. 

2 Background 

2.1 Notation 

The following notation is used in this paper. Let L(C n ,C m ) be the set of all linear operators 
from C n to C m and let L(C n ) := L(C n ,C n ). Next, let Pos(C n ) C L(C n ) be the set of all positive 
semidefinite operators on C n . Let ||Af|| max : = maxjj \Mij\ denote the largest entry of M G L(C n ) 
in absolute value. Finally, for any natural number re, let [re] := {1, . . . , n} and let I n be the re x n 
identity matrix. 

2.2 Separable and LOCC measurements 

A A:-outcome POVM measurement (or simply a measurement) on an re-dimensional state space is 
a set of operators {E%, . . . , E^} C Pos(C n ) such that Yli=i = In- The operators Ei are called 
POVM elements. The probability of obtaining outcome i upon measuring state p is Tr (Eip). 

When it is necessary to keep track of the post-measurement state, it is more convenient to use a 
non- destructive measurement. Such a measurement is specified by a set of measurement operators 
{M±, . . . , .Mfc} C L(C n , C m ) for some finite m where YLi=i Mi = I n . The probability of obtaining 
outcome i upon measuring state p is Tr(Mj Mip) and the rre-dimensional post-measurement state 
is MipMj /Tr (MjMip). 

Note that a non-destructive measurement followed by discarding the post-measurement state 
corresponds to a POVM measurement with elements Ei = MjMi. 
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2.2.1 Separable measurements 

Definition 1. A measurement £ = {E\, . . . , E k } on a bipartite state space C dA (g) C dB is separable 
if all POVM elements E^ are separable, i.e., 

Ei = Yl E f® E f C 1 ) 

3 

for some Ef £ Vos{C dA ) and Ef £ Pos(C ds ). 

Note that the above definition is equivalent to saying that M is obtained from a measurement 
with product POVM elements, followed by classical post-processing (coarse graining). 

2.2.2 LOCC measurements 

Informally, a bipartite n-outcome LOCC measurement £ consists of the two parties taking finitely 
many turns (called rounds) of applying adaptive non-destructive measurements to their state spaces 
and exchanging classical messages. This is followed by coarse graining all measurement records into 
n bins, each corresponding to one of the n outcomes of £. 

Let us describe such a protocol £ more formally, adopting notation similar to that of [BDF + 99]. 
Let A denote the empty string, corresponding to no message being sent. The protocol begins when 
one of the parties, say Alice, applies a non-destructive measurement 

A(A) = {A 1 (A),...,A k{A) (A)} (2) 

to her state space and communicates the round 1 measurement outcome mi £ [&(A)] to Bob. Then, 
depending on the value of mi received, Bob applies a non-destructive measurement 

B(mx) = {BiK), . . . ,B k{mi) {mi)} (3) 

to his state space and communicates the round 2 measurement outcome mi £ [fc(mi)] to Alice. 
The protocol proceeds with the two parties taking finitely many alternating turns of a similar form, 
where the non-destructive measurement applied at round t depends on the measurement record 
m = (mi, . . . , mt_i) accumulated during the previous rounds. 

Let m be the measurement record after the execution of the first t rounds of the protocol. Then 
the measurement operator that Alice and Bob have effectively implemented is a product operator 
A m <8> B m , where 1 

A m := A mt ^ 1 (m 1 , . . . ,m t _ 2 ) • . . A mz (mi, m 2 )A mi (A), (4) 
B m := B mt (mi, . . .,m t -i) ■ ■ ■ B nH (mi,m 2 ,m 3 )B m2 (mi). (5) 

Alice and Bob may choose to terminate the protocol depending on the measurement record 
obtained. At this point they must output one of the n outcomes of the LOCC measurement £ that 
they are implementing. If L{k) is the set of all terminating measurement records corresponding to 
outcome k £ [n], then the fcth POVM element of £ is given by 

E k := A l A m ® B m B m . (6) 

m£L(k) 

Since each E k is separable, any LOCC measurement is separable. 

1 Here we assume for simplicity that t is even; in the odd case the operators Am and B m can be defined similarly. 
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2.2.3 Finite and asymptotic LOCC 



We consider two scenarios: when a measurement can be performed in a finite number of rounds or 
asymptotically. 

Definition 2. We say that a measurement £ can be implemented by (finite) LOCC if there exists a 
finite-round LOCC protocol that, for any input state, produces the same distribution of measurement 
outcomes as £. 

Definition 3. We say that a measurement £ can be implemented by asymptotic LOCC if there 
exists a sequence V\,Vi, ■ ■ ■ of finite-round LOCC protocols whose output distributions converge to 
that of £. 

The exact implementation scenario is not practical since any real-world device is susceptible 
to errors due to imperfections in implementation. However, proving that a certain task cannot be 
performed asymptotically is considerably harder than showing that it cannot be done (exactly) by 
any finite LOCC protocol. 

2.2.4 LOCC protocol as a tree 





Bi(l) 



^i(A) A 2 (A) A 3 (A) 

(7) m 

A„( I 

5 2 (1) 




A 2 (A)®I 



A 3 (A)(8)J 




1,1 



Ai(A) ® #i(l) 




Ai(l,2) A 2 (l,2) 



1,2,1 




A 1 (l,2)A 1 (A)®B 2 (l) 



A 2 (l,2)Ai(A)®B 2 (l) 



Figure 1: Tree structure of a three-outcome LOCC measurement. In round one Alice performs a three- 
outcome non-destructive measurement .4(A); in round two, upon receiving message "1", Bob performs a 
two-outcome non-destructive measurement 6(1) and upon receiving message "2" or "3" he terminates the 
protocol; in round three, upon receiving message "1" , Alice terminates the protocol and, upon receiving 
message "2", she performs a two-outcome non-destructive measurement ^4(1,2). All nodes are labeled by 
the accumulated measurement record. The corresponding measurement operator is given below each leaf. 



We represent an LOCC measurement protocol as a tree (see Figure 1). The protocol begins 
at the root and proceeds downward along the edges. Each edge represents a certain measurement 
outcome obtained at its parent node, and leaves are the nodes where the protocol terminates. 
The set of all leaves is partitioned into subsets, each corresponding to an outcome of the LOCC 
measurement being implemented. 
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A path from the root to a leaf is called a branch. There is a one-to-one correspondence between 
the branches and the possible courses of execution of the LOCC protocol. Likewise, there is a one- 
to-one correspondence between the nodes of the tree and the accumulated measurement records. 

The measurement at node u is the measurement performed by the acting party once the protocol 
has reached node u. In contrast, the measurement operator corresponding to node u is the measure- 
ment operator that has been implemented upon reaching node u. For example, consider the node 
(1, 2). The measurement at node (1, 2) is given by the POVM {Ax(l, 2), ^(L 2)}, whereas the mea- 
surement operator corresponding to the node (1, 2) is given by ^4i(A) (g>£>2(l)- As another example, 
the measurement operators corresponding to the leaves are exactly the measurement operators of 
the LOCC protocol prior to coarse graining. 

2.3 Bipartite state discrimination problem 

The goal of this paper is to investigate the limitations of two-party LOCC protocols for the task 
of bipartite quantum state discrimination, which is as follows: 

Let S = {\ipi), • • • j iV'n)} C C dA C B be a known set of quantum states. Suppose that k € [n] is 
selected uniformly at random and Alice and Bob are given the corresponding parts of state \ipk) £ S. 
Their task is to determine the index k by performing a measurement on this state. 

A case of special interest is when S is an orthonormal product basis, i.e., each = \cti)\{3i) 
for some orthonormal bases \ai) £ C dA and € C dB . Such states can be perfectly discriminated 
by a separable measurement £ with POVM elements 

E i :=\a i ){a i \®\p i ){p i \. (7) 

However, this measurement cannot always be implemented by finite [WH02, GV01] or even asymp- 
totic LOCC [BDF + 99]. In such cases we say that S possesses nonlocality (without entanglement). 

2.4 Previous results 

The first example of an orthonormal product basis of bipartite quantum states that cannot be 
perfectly discriminated by (even asymptotic) LOCC was given in [BDF + 99]. This is a striking 
illustration of the difference between the power of LOCC and separable operations. Furthermore, 
[BDF + 99] quantifies the information deficit of any LOCC protocol for discriminating these states. 
This result has been a starting point for many other studies on state discrimination by LOCC, 
with the ultimate goal of understanding LOCC operations and how they differ from separable ones. 
We briefly describe some of the directions that have been explored. Unless otherwise stated, these 
results refer to the discrimination of pure states with finite LOCC. 

First consider the problem of discriminating two states without any restrictions on their dimen- 
sion. Surprisingly, any two orthogonal (possibly entangled) pure states can be perfectly discrimi- 
nated by LOCC, even when they are held by more than two parties [WSHVOO]. Furthermore, opti- 
mal discrimination of any two multipartite pure states can be achieved with LOCC both in the sense 
of minimum error probability [VSPM01] and unambiguous discrimination [CY01, CY02, JCY05]. 
Recently this has been generalized to implementing an arbitrary POVM by LOCC in any 2- 
dimensional subspace [Crol2]. 

Many authors have considered the problem of perfect state discrimination by finite LOCC. In 
particular, the case where one party holds a small-dimensional system is well understood. Reference 
[WH02] characterizes when a set of orthogonal (possibly entangled) states in C 2 (g> C 2 can be 
perfectly discriminated by LOCC. A similar characterization for sets of orthogonal product states 
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Ill C 3 <g> C 3 has been given by [FS09]. In addition, [WH02] characterizes when a set of orthogonal 
states in C 2 (g> C n can be perfectly discriminated by LOCC when Alice performs the first nontrivial 
measurement. It is also known that ^-rotated domino states cannot be perfectly discriminated by 
LOCC (unless 6 = 0) [GV01]. Furthermore, the original domino states have inspired a construction 
of n-partite d-dimensional product bases that cannot be perfectly discriminated with LOCC [NC06]. 

The role of entanglement in perfect state discrimination by finite LOCC has also been consid- 
ered. It is not possible to perfectly discriminate more than two Bell states by LOCC [GKR + 01]. 
In fact, the same is true for any set of more than n maximally entangled states in C n <8> C n [Nat05]. 
Multipartite states from an orthonormal basis can be perfectly discriminated by LOCC only if it 
is a product basis [HSSH03]. Also, no basis of the subspace orthogonal to a state with orthogonal 
Schmidt number 3 or greater can be perfectly discriminated by LOCC [DFXY09]. On the other 
hand, any three orthogonal maximally entangled states in C 3 (g) C 3 can be perfectly discriminated 
by LOCC [Nat05]. In fact, if the number of dimensions is not restricted, one can find arbitrar- 
ily large sets of orthogonal maximally entangled states that can be perfectly discriminated by 
LOCC [Fan04]. Contrary to intuition, states with more entanglement can sometimes be discrim- 
inated perfectly with LOCC while their less entangled counterparts cannot [HSSH03]. Generally, 
however, a set of orthogonal multipartite states S C C D can be perfectly discriminated with LOCC 
only if \S\ < ^^|y> where d(S) measures the average entanglement of the states in S [HMM + 06]. 

It is known that local projective measurements are sufficient to discriminate states from an 
orthonormal product basis with LOCC [DR04, CL04]. Moreover, there is a polynomial-time (cubic 
in maxjc^, ds}) algorithm for deciding if states from a given orthonormal product basis of C dA 
C dB can be perfectly discriminated with LOCC [DR04]. The state discrimination problem for 
incomplete orthonormal sets (i.e., orthonormal sets of states that do not span the entire space) seems 
to be harder to analyze. However, unextendible product bases might be an exception (although 
commonly referred to as "bases" these are in fact incomplete orthonormal sets). It is known 
that states from an unextendible product basis cannot be perfectly discriminated by finite LOCC 
[BDM + 99]. In fact, the same holds for any basis of a subspace spanned by an unextendible product 
basis in C 2 <%> C 2 <g> C 2 [DXY10] . Curiously, there are only two families of unextendible product 
bases in C 3 ® C 3 , one of which is closely related to the domino states [DMS + 03]. 

The problem of state discrimination with asymptotic LOCC has been studied less. It is known 
that states from an unextendible orthonormal product set cannot be perfectly discriminated with 
LOCC even asymptotically [DR04]. Reference [KKB11] gives a necessary condition for perfect 
asymptotic LOCC discrimination, and also shows that for perfectly discriminating states from an 
orthonormal product basis, asymptotic LOCC gives no advantage over finite LOCC. The latter 
result implies that the algorithm from [DR04] also covers the asymptotic case. On the other hand, 
even in some very basic instances of state discrimination it remains unclear whether asymptotic 
LOCC is superior to finite LOCC (see [DFXY09, KKB11] for specific sets of states). 

Another line of study originating from [BDF + 99] aims at understanding the difference be- 
tween the classes of separable and LOCC operations. To this end, [Cohll] constructs an r-round 
LOCC protocol implementing an arbitrary separable measurement whenever such a protocol ex- 
ists. A different approach is to exhibit quantitative gaps between the two classes. To the best 
of our knowledge, only two quantitative gaps other than that of [BDF + 99] are known. References 
[KTYI07, Koa09] demonstrate a gap between the success probabilities achievable by bipartite sepa- 
rable and LOCC operations for unambiguously discriminating 1 00) from a fixed rank-2 mixed state. 
The largest known difference between the two classes is a gap of 0.125 between the achievable suc- 
cess probabilities for tripartite EPR pair distillation [CCL11]. Moreover, as the number of parties 
grows, the gap approaches 0.37 [CCL11]. 

At a first glance one might think that the nonlocality without entanglement phenomenon is 
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related to quantum discord. However, the quantum discord value cannot be used to determine 
whether states from a given ensemble can be discriminated with LOCC [BT10]. 

Finally, if a set of orthogonal (product or entangled) states cannot be perfectly discriminated 
by LOCC, one can measure their nonlocality by considering how much entanglement is needed to 
achieve perfect discrimination [C0I1O8, BBKW09]. 

3 Framework 

In this section we introduce a framework for proving lower bounds on the error probability of any 
LOCC measurement for discriminating bipartite states from a given set 

S:={|Vi},...,|Vn>} cC dA ®C dB . (8) 

We make no assumptions about the states In particular, they need not be product states or 
be mutually orthogonal. 

From now on, V denotes an arbitrary LOCC protocol for discriminating states from S. In rough 
outline our argument proceeds as follows: 

1. We modify V so that it can be stopped when a specific amount of information e has been 
obtained (see Section 3.1). This is done by terminating the protocol prematurely and possibly 
making the last measurement less informative (see Section 3.2). 

2. When the information gain is e, we lower bound a measure of disturbance (defined in Sec- 
tion 3.3) by rje for some constant r/ (see Section 3.4). 

3. We show that at least two of the possible initial states have become nonorthogonal at this 
stage of the protocol, and we infer a lower bound on the error probability of V (see Section 3.5). 

Our framework reuses some ideas of the original approach [BDF + 99]. However, instead of mu- 
tual information, we quantify how much an LOCC protocol has learned about the state using error 
probability. This allows us to replace the long mutual information analysis in the original paper 
with a simple application of Helstrom's bound. The idea of relating information gain and distur- 
bance also comes from [BDF + 99]. Here, we analyze this tradeoff using the nonlocality constant 
(see Definition 3.4) which can be applied to any set of states. In Section 4 we give a method for 
lower bounding the nonlocality constant that applies specifically when S is an orthonormal basis of 
C dA <g>C ds . In Section 5 we apply this method for the domino states and some other related bases. 

3.1 Interpolated LOCC protocol 

Consider an arbitrary node in the tree representing the protocol V . Let m be the corresponding 
measurement record and let A <g> B denote the Kraus operator that is applied to the initial state 
when this node is reached. Note that the output dimensions of operators A and B could be 
arbitrary. 

The initial state \ipk) yields measurement record m with probability 

p(m|^ fe ):=Tir[(A®B)t(A®B)|^)(^ fc |] = {if> k \(a ® 6)|Vft> (9) 

where a := A* A E Pos(C dA ) and b := B^B € Pos(C ds ). Note that we need not concern ourselves 
with the arbitrary output dimensions of A and B from this point onward. We use Bayes's rule and 
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the uniformity of the probabilities p{ipk) to obtain the probability that the initial state was \ipk) 
conditioned on the measurement record being m: 



p(ipk\m) 



p(4>k)p(m\il; k ) 



{ip k \{a®b)\i) k ) 



ELiPW'jMHV'j) Ej=i(V>jl(a ® b)\ipj 



(10) 



At the root, the measurement record m is the empty string and p(ipk \ m ) = ^ f° r au k. As we proceed 
toward the leaves, these probabilities fluctuate away from -. For example, if V discriminates the 
states perfectly, the distribution reaches a Kronecker delta function. 
For a given node m let us define 



; (m) := m&xp(ipk\ m )- 
fceFnl 



(11) 



Let e := 



p max (m) — Then e characterizes the uniformity of the distribution and thus 

the amount of information learned about the input state. The next theorem shows that we can 
modify the protocol V so that it can be stopped when some but not too much information has been 
learned. While this idea originates from [BDF + 99], we use a specific result from [KKB11]. 

Theorem 1 (Kleinmann, Kampermann, BruB [KKB11]). LetV be an LOCC protocol for discrim- 
inating states from a set S of size n. For any e > there exists an LOCC protocol V e that has the 
same success probability as V, but each branch ofV £ has a node m such that either 



1 



Pmax(m) = - +e 



1 



or 



n 



Pmax(m) < — he and m is a leaf ofV. 



(12) 



Proof idea. Let u be a node in the protocol tree of V and let V\ , . . . , v m be the children of u. Assume 
that for some i we have 



Pmax(-u) < - + S < PmaxK), 
n 



(13) 



which means that the measurement outcome corresponding to the edge (n, Vi) is too informative. 
To rectify this, we break up the measurement at node u into two steps. We represent the outcomes 
of the first measurement by new nodes v\, . . . ,v m while the outcomes of the second measurement 
lead to the original nodes vi,...,v m (see Figure 2). 



Protocol V: 




Protocol V E : 




Figure 2: The protocol tree before (left) and after (right) splitting the measurement at node u into two 
steps. (The graph on the right has been condensed for clarity, but it can be expanded into a tree by making 
a new copy of subtree Ti for each incoming arc in Vi.) The amount of information learned in the first step 
is controlled by diluting the measurement operators, and the purpose of the second step is to complete the 
original measurement. The dotted line corresponds to the end of stage I (see Definition 4). 
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The first measurement interpolates between a completely uninformative trivial measurement 
and the original measurement at u. The interpolation parameters are chosen so that p max (vi) = 
for all i that satisfy Equation (13). The second measurement depends on the outcome of the first 
measurement. It produces the same set of post-measurement states as the original measurement at 
u. Moreover, the total probability of obtaining each state is the same as in the case of the original 
measurement. After this we proceed according to the original protocol. 

Protocol V e is obtained from V by considering all branches of V and performing the above 
procedure at the closest node to the root that has a child satisfying Equation (13). For more 
details see [KKB11]. □ 

In the context of state discrimination, the possibility of interpolating a protocol to obtain some 
but not too much information is what distinguishes LOCC measurements from separable ones. In 
particular, a separable measurement for a set of states that cannot be distinguished by asymptotic 
LOCC cannot be divided into two steps, with the first yielding information precisely e and the 
second completing the measurement (further details will be provided in a manuscript currently in 
preparation). 

3.2 Stopping condition 

To control how much information the protocol has learned, we fix some e > and stop the execution 
of V £ when we reach a node m that satisfies the conditions in Equation (12). 

Definition 4. We say that stage I of the protocol V £ is complete at the earliest point when Equa- 
tion (12) is satisfied. 

We choose e < n ^_^ in our analysis. Operationally, this means that none of the n states has 
been eliminated at the end of stage I, since 

min p(%[)k\m) > 1 — (n — l)p max (m) > (n — l)e > 0. (14) 

ke[n] n 

This allows us to use Helstrom's bound to lower bound the probability of error (see Section 3.5). 
It also ensures that the disturbance measure 5s(a <S> b) introduced in Section 3.3 is well defined at 
m. All constraints imposed on the distribution p(ipk\ m ) are summarized in Figure 3. 

p(ipk\m) 



Mil 



(n - l)e 



Figure 3: Probability distribution p(ij)k\rn) at the end of stage I. For all k we have — + e > p(ipk\m) > 
(n — l)e > where the first inequality is tight for some fc. 

Since the error probability of the protocol V e is a weighted average of error probabilities of 
individual branches, it suffices to lower bound these individual error probabilities. For any branch 
that terminates without a node satisfying 

Pmax(m) = - +£, (15) 
n 
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we can put a large lower bound on the error probability. In particular, for the optimal choice 



2 1 



3 n(n— 1) 



of Theorem 2 with n > 2, 

/1\ 121 1 

Perror(m-) > 1 -p max m >1--+£=1 -— — > -, (16) 

\ra / n 3 n(n — 1) 6 

which is much higher than the lower bound we obtain for other branches. We now consider the 
remaining case where stage I ends with a node satisfying Equation (15). 

3.3 Measure of disturbance 

Now we show that at least two possible post-measurement states (A (g) B)\tpi) and (A (8> B)\ipj) 
are nonorthogonal at the end of stage I, and lower bound their overlap quantitatively. Assuming 
that the initial state was E S, the normalized post-measurement state at the node with 
corresponding measurement operator A <g> B is 

where a := A^A and b '■= B^B. Note that I ( a ® b)\ipi) > for all i E [n] because, from 
Equations (14) and (10), < min fce[n ] p(ip k \m) = min fce[n] ^^(^a^t)^) • 

Definition 5. The disturbance caused by the operator a®b on the set of states S is defined as 
5 s (a®b) := max max , Wi)! (18) 

Note that 5s(a®b) measures the nonorthogonality of the post-measurement states. If the initial 
states were orthogonal then 5 s (a ® 6) indeed characterizes the disturbance caused by a ® b. 

Since (<f>i\<f>j) can be expressed in terms of the operators a = A^A and b = B^B, from now on 
we no longer explicitly use the measurement operators A and B. 

3.4 Disturbance/information gain trade-off 

Now we define the nonlocality constant and show that it relates 5 (the disturbance caused at the 
end of stage I) to e (the amount of information learned). 

Definition 6. The nonlocality constant of S is the supremum over all rj such that for all a E 
Pos(C dA ), b E Pos(C ds ) and for all i satisfying {ipi\(a (g> b)\ipi) ^ 0, 

f max ke[n] (ip k \(a ® b)\ip k ) 1\ 

^ " — / , i/ ~r\m < 5s{a ® 6) . (19) 

V E i6 [„]<^l(o®b)l^> «/ 

Equivalently, if Gy := | (a (g) fr)!^') f° r E [n] then 

' |G, 



r/ := inf < 

a,b 



l 3 

maxj^j 



y/ GaG 



33 



max fc G fefc 1 



(20) 



33 



n 



where the infimum is over all a E Pos(C ciA ) and b E Pos(C rfs ) such that Gu ^ for all i G [n]. 

Recall from Section 3.2 that we stop the LOCC protocol at the end of stage I in a node m 
where the condition in Equation (15) is satisfied for some e E (0, n ^-i) )~ °® ^ be the operator 
corresponding to node m and let 5 := 5s (a (g) b) be the disturbance caused. 
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Lemma 1 (Disturbance/information gain trade-off). The amount of information e learned at the 
end of stage I lower bounds the disturbance 5 as 

r]e<5 (21) 

where r\ is the nonlocality constant of S (see Definition 3.4)- 

Proof. This immediately follows from the definitions of e and n: 

( ,. , v 1\ f m&x ke[n] (ij k \(a ® b)\ip k ) 1\ 
7]E = r) maxp \m) } = V\ i , r, - , r, , \ ]<o (22) 

where we have used Equations (15), (10), and (19). □ 



3.5 Lower bounding the error probability 

In this section we use Lemma 1 to lower bound the error probability of any LOCC measurement 
for discriminating states from the set S. 

Note that Equation (21) together with the definition of 5 implies that at the end of stage I there 
are two distinct post-measurement states \4>i) and \<j)j) such that 

=d> V e. (23) 

As discussed in Section 3.2, our choice of e guarantees that p(ipi\m) and p(ifjj\m) are both strictly 
positive. Thus we can use the following result to lower bound the error probability: 

Fact (Helstrom bound [Hel76, pp.113]). Suppose we are given state \<&q) with probability qo and 
state |<I>i) with probability q% = 1 — q$. Any measurement trying to discriminate the two cases errs 
with probability at least 

Q(qo,qi,S) :=l{l-V^-^oqiS 2 )>q qiS 2 , (24) 

where 5 = |($o|^l}| is the overlap between the two states, and the inequality follows from 1 — 
Vl - x 2 > \x 2 for x £ [0, 1]. 

As e increases, the disturbance (thus the overlap between some and \4>j)) increases, but the 
lower bound on the probabilities p{ipi\m) and p{il)j\m) decreases. The choice e = f ^fazrn gi yes a 
lower bound on the error probability as follows. 

Theorem 2. Let S be a set of quantum states in <C dA ®C dB of size n > 2. Any LOCC measurement 
for discriminating states drawn uniformly from S errs with probability 

2 r? 2 

Pcrror > ^ ^5 (25) 

where rj is the nonlocality constant of S (see Definition 3.4)- 

Proof. At the end of stage I there are two post-measurement states l^o) an d l^i) with overlap 5. 
Let po and p\ be the posterior probabilities of these states. To lower bound the error probability 
of V e (thus that of V), we give Alice and Bob extra power at this point: 

• if the actual input state does not lead to |<I>o) or we assume that Alice and Bob succeed 
with certainty; 
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• otherwise Alice and Bob are allowed to perform the best joint measurement to discriminate 
the states |3>o) an d 

For fixed e and probabilities po and pi, we can lower bound the error probability by the following 
expression: 

P(P0,Pi,e):=(p + p 1 )-Q(^, ¥ ^,8). (26) 
Using Equation (24) and the inequality 5 > r/e from Lemma 1, we get that 

P( Po ,Pi,e)>^-(ne) 2 . (27) 
Po+Pi 

Recall that we stop the protocol at a point where we are guaranteed that < e < n ^-i) an d> 
by Equations (14) and (15), 

--(n-l)s< Pi <-+E (28) 
n n 

for all i. Given these constraints on po and p\, we can choose the e that maximizes P(po,pi, e) and 
guarantee that the error probability in the branch of the LOCC protocol being considered satisfies 

Pcrror > max min P(p ,pi,e). (29) 

£e (°.^T)) Po,Pi6[i-(n-l)ei+ £ ] 

From Equation (27) we get 

Pcrror > niax mm (r]£) . (SO) 

^(0,^) POlPie [l-(„-l) £ l +£ ] Po+Pi 

The minimum is attained when po = p\ = ^ — (n — l)e (i.e., the probabilities are equal and as 
small as possible), so the problem simplifies to 

1/1, N \, , 2 2 rj 2 2t) 2 

Perror > max - - - (n - 1 £ J fa e) >^~37 TT2>^;— (31) 

\ ' n(n — 1) / 

where the value 

2 1 

e = ; r 32 

3n(n-l) V 7 

achieves the maximum. □ 

Theorem 2 shows that any LOCC protocol for discriminating states from S errs with probability 
proportional to rf, justifying the name "nonlocality constant." 



4 Bounding the nonlocality constant 

The framework described in Section 3 reduces the problem of bounding the error probability for 
discriminating bipartite states by LOCC to the one of bounding the nonlocality constant r] (see 
Theorem 2). This reduction holds for any set of pure states S. In this section we assume that S 
is an orthonormal basis of C dA <S> C rfs and provide tools for bounding the nonlocality constant. In 
particular, we bound r] in terms of another quantity that we call "rigidity" . 

For the remainder of the paper we represent pure states from C dA <8> C dB using "tiles" in a 
cLa x d>B grid. We first introduce some notations related to tilings in Section 4.1. Then we define 
rigidity and relate it to the nonlocality constant r\ in Section 4.2. Section 4.3 provides a tool, the 
"pair of tiles" lemma, that we use to bound rigidity for specific sets of states in Section 5. 



14 



4.1 Definitions 



Given a fixed orthonormal basis {\i) : i G [d]}, define the support of a pure state |^>) E C d as 



If |V>) G C dA ® C ds then supp C [d^] x [rffl]- Consider [oIa] x [d#] as a rectangular grid of size 
oIa x cIb- Any region that corresponds to a submatrix of this grid is called a tile. More formally, a 
i«Ze is a subset T C [d^] x [d#] such that T = R x C for some i? C [d^] and C C [dfl]. (Note that 
a tile is not necessarily a contiguous region of the grid.) We use rows(T) = R and cols(T) = C to 
denote the rows and columns of this tile, respectively, and we use \T\ to denote the size or the area 
of T. If |^>) = | a) |/3) is a product state, then supp \ ip) = supp |a) x supp \ j3) and thus supp \ip) is a 
tile, which we call the tile induced by 

We say that an orthonormal set of product states S C C dA ® C dB induces a tiling of a d^ x rf B 
grid if the tiles induced by the states in S are either disjoint or identical. Note that if S is an 
orthonormal basis of C dA <S> C dB , then a tile of area L is induced by L states that form a basis of 
that tile. In a domino-type tiling, every tile has area 1 or 2. 

For a given tiling T of a d^ x d# grid let us define the corresponding row graph as follows: its 
vertex set is [d^] with two vertices i and j adjacent if and only if there exists a column c such that 
(i, c) and (j, c) belong to the same tile. The column graph of a tiling is defined similarly. We say 
that a tiling is irreducible if its row graph and its column graph are both connected. The diameter 
of the tiling T is the maximum of the diameters of its row and column graphs. See Figure 4 for an 
example. 



Figure 4: A domino-type tiling and the corresponding row and column graphs. This tiling is irreducible and 
has diameter two. 

Without loss of generality we consider only irreducible tilings. Reducible tilings can be broken 
down into several smaller components without disturbing the underlying states. To do this, both 
parties simply perform a projective measurement with respect to the subspaces corresponding to 
the different components of the row and column graphs. 

Note that in general, a tiling is not invariant under local unitaries. In particular, the irre- 
ducibility of the tiling induced by a given set of states is a basis-dependent property. The most 
extreme example of this phenomenon is the case of the standard basis. It induces a completely 
reducible tiling that consists only of 1 x 1 tiles. However, if both parties apply a generic local 
unitary transformation, the resulting tiling consists only of a single tile of maximal size. 

4.2 Lower bounding the nonlocality constant using rigidity 

In this section we assume that S is an orthonormal basis of C dA (g>C ds (so in particular, n = dj^ds) 
and discuss a particular strategy for lower bounding r\ for such S. We apply this strategy to several 
sets of orthonormal product bases in Section 5. 



supp |V) := {i G [d] : (#) ^ 0} . 



(33) 
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We bound rj (quantifying a disturbance/strength tradeoff) by considering a quantitative prop- 
erty of the set S called rigidity. Intuitively, we call a measurement operator strong if it is far from 
being proportional to the identity matrix; a set of states S is rigid if there exists a strong measure- 
ment that leaves the set undisturbed. We formalize this as follows (recall that ||-|| max denotes the 
largest entry of a matrix in absolute value): 

Definition 7. For an orthonormal basis S, if there is a constant c such that for all a G Pos(C rfA ), 
b G Pos(C ds ) and for all i such that {ipi\(a ® b)\ipi) / 0, 



a® b I 



< c-5 s (a®b), (34) 



Tr(a (g) b) n 

we say S is c-rigid, or c is an upper bound on the rigidity of S. 

When S is rigid, the states can remain unchanged despite application of a strong measurement. 
For example, a tensor product basis is not c-rigid for any finite c (i.e., such a basis is arbitrarily 
rigid). In contrast, if c is small, then any strong measurement disturbs the set S, and Equation (34) 
quantifies how weak a measurement operator a <%> 6 must be for the disturbance 5s (a (g) b) to be 
small. 

We now relate upper bounds on the rigidity of S to lower bounds on its nonlocality constant: 
Lemma 2. Let S be an orthonormal basis of C dA (g) C dB . If S is c-rigid then 

i £ 7i- (35> 

where L is the size of the largest tile corresponding to states in S. 

Proof. If S is c-rig id, then for any a G Pos(C dA ) and b G Pos(C ds ) (such that (ip k \(a ® b)\ip k ) + 
for all k G [n]), we have 

°®J^ -- = cM-* s ( ®6) (36) 
Ir(a (8) b) n 

for some Hermitian matrix M G L(C dA (8> C dB ) with ||-/Vf|| max < 1. From this we get 

max(V>fc| "® J\ J V'fc) - - = cmax(^|M|V>fc) • 5 s (a ® b) (37) 
ke[n] ir{a (£> b) n fce[n] 

<cL-5 s (a®b). (38) 

By the definition of r/ (Equation (19)) and the fact that Tr(a (8) b) = X]je[n](^jl( a ® b)\ipj) for any 
orthonormal basis S, we get the desired inequality. □ 

Putting Lemma 2 and Theorem 2 together gives the following: 

Theorem 3. Let S be an orthonormal basis of C dA C dB . If S is c-rigid then any LOCC mea- 
surement for discriminating states from S errs with probability 

2 1 

Pemi ~ 27 (cL)%5 ( 39 ) 
where L is the size of the largest tile of S. 
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4.3 The "pair of tiles" lemma 

In this section we present a lemma that serves as our main tool for bounding rigidity. 

= U\i) for i G [m] and \ifij) := V\j) for 



Lemma 3. Let U G \J(m),V G U(n), and define \<pi 
j G [n]. Then for any M G L(C n ,C m ) we have 



/ mn ■ max\(ipi\M\'iljj)\ >max|Mjy|. 
i,j k,l 



(40) 



The main idea of the proof is that a unitary change of basis can only increase the largest entry 
of a vector by a multiplicative factor depending on the dimension of the vector. 



Proof. Let us define a mapping vec: L(C n , C m ) — > C n (g> C' m as 

vec: H- \i)\j) 



(41) 



for i G [m] and j G [n] and extend it by linearity over C. One can check that vec(AX B) = 
(A <g> B T ) vec(X). Using this and basic inequalities between the 2-norm and the oo-norm, we get 



max|(^|M|^)| 
i,3 



yec{J2(<Pi\M\^)\i)(j\ 
vec(j2(i\U^MV\j)\i)(j\ 



vec 



> 



> 



'run ' 
1 . 
Imn ' 
1 | 
'mn 
1 



MIL 

(C/t ® y T )vec(M)|| 2 



® V T )vec 

1 



|vec(M)|| 2 
I vec (M) II 



max | Mj 
yjmn k,i 

as desired. 

Let us restate Lemma 3 using the language of tilings: 



kl 



(42) 

(43) 

(44) 
(45) 

(46) 
(47) 
(48) 
(49) 
□ 



Lemma 4. Let R\,R2 C [o!^] x [fig] 6e too arbitrary regions of a dA x &b grid, and {|yi)}i=i and 
{iV'i)}^! C C dA (giC ' 3 6e £/ieir 6ases (Tiere \(fi) and \ipj) need not be product states). Then for any 
matrices a G L(C dA ) and b G L(C rfs ) we have 



y/\Ri\-\R,2\ max\{<pi\(a ® b)\i})j)\ > max \a r 
id ' ' (ri,ci)ei?! 



J C\C2 I 



(50) 



This follows from Lemma 3 by restricting |^j) and to regions R\ and i?2, respectively, and 
choosing M to be a submatrix of a ® b with rows determined by R\ and columns by ife- 
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Proof. For t G {1, 2} let us enumerate the cells of region Rt by integers from {1, . . . , \Rt\} arbitrarily, 
and let (rj(i), be the coordinates of the ith. cell of region Rt- Let 

\Rt\ 

n t :=J2\i)(n(i),ct(i)\ (51) 

i=l 

be a linear operator that restricts the space C dA <X> C dB to region R t . Then |<^) := IIi|y>j) is 

the restriction of \(fi) to region R\ and |^) := II2 1 "04 ) is the restriction of \tfji) to R2. Also, let 
M := IIi(a® £>)l4. 

Note that for all i G {1, . . . , |i?i|} we have njTfi|c^j) = \<pi) since the support of \(fi) lies 

entirely within region R\ and njlli is the projection onto R±. Similarly, 112112 1 V^j) = f° r an 
j G {1,...,|# 2 |}. Hence 

(y>i|(a®6)|^) = (^|n t i n 1 (a»6)4n 2 |^) = (<Pi\M\^) (52) 
for all i and j. Finally, we apply Lemma 3 to (1^)}^', {\tp'j)}^l, and M: 

#i| • l^j max I ( Vi |(a® 6)|^)| = ylWl max|(^|M|#)| 

> max\M k i\ 
k.l 



max 

k,l 



|(n(fc)|a|r 2 (0)|-|<ci(fc)|6b(0>| 



max |(fc|IIi(a(g> &)rLj, 



— max |tt rir2 | • |6 ClC2 | 
(n,ci)eiii 

(»-2,C2)ei?2 

and the result follows. □ 

When regions i?i and R2 are two distinct tiles from the tiling induced by S, we can use Lemma 4 
to get the following result: 

Lemma 5 ("Pair of tiles" Lemma). Let T\ and T 2 be two distinct tiles in the tiling induced by S, 
and let a G Pos(C dA ) and b G Pos(C ds ). Then 

y/\Ti\ • \T 2 \ 5 s (a ® b) Tr(a ® b) > \ I (53) 

for any rt G rows(Tt) and Ct G cols(Tt) where t G {1,2}. 

Proof. We relax the inequality in Lemma 4 by observing that 

max^j I (^1(0(8)6)1^) I > max^- |(^,|(a® b)|V>j)| , , 

5 l|a®&lloo ~ Tr(a®6) 1 ' 

which easily follows from the definition of d~s(a ® b) in Equation (18). □ 

Note that the tiles T\ and T 2 in Lemma 5 have to be distinct since the maximization in the 
definition of <5s(a<X> b) is performed only over pairs of distinct states. This lemma will be used later 
to bound the off-diagonal entries of a ® b (see Figure 5). 
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r-2 



ci 



+ 



C-2 



Cl = c 2 

ri 

r2 - 



Figure 5: Whenever (ri, ci) and (?~2, C2) belong to different tiles (left), Lemma 5 can be used to upper bound 



the off-diagonal entry a rir2 ■ i 
result cannot be applied directly. 



of a £§> b. When both coordinates correspond to the same tile (right), this 



5 Domino states 

In this section we use the framework introduced earlier to give a lower bound on the error probability 
of any LOCC measurement for discriminating states from certain bipartite orthonormal product 
bases known as domino states. This provides an alternative proof of the quantitative separation 
between LOCC and separable measurements first given in [BDF + 99] as well as generalizations to 
states corresponding to other domino-type tilings and a rotated version of the original domino 
states. 

Bob 

|0> |1) |2) 

|0) 

Alice |1> 
|2> 

Figure 6: The tiling induced by states from Equations (55-59). 



2,3 



6,7 



8,9 



4,5 



5.1 Definition 

The following orthonormal product basis is known as the domino states: 



h/h) = |l>|l>, (55) 

hfe) = |0)|0 + l>, |^ 3 ) = |0)|0-1), (56) 

|V>4> = |2)|l + 2>, |V>5> = |2>|l-2>, (57) 

|V>6> = |1 + 2)|0), |V>7> = |1-2)|0), (58) 

|V> 8 > = |0 + 1)|2), ^9) = |0-1)|2), (59) 



where \i±j) := (\i)±\j))/>/2. In [BDF+99] it was shown that any LOCC protocol for discriminating 
these states has information deficit at least 5.31 x 10 -6 (out of log 2 9 ~ 3.17) bits. 

In [BDF + 99] the authors also consider a family of orthonormal product bases, the so-called 
rotated domino states, which are parametrized by four angles < 01,82,63,64 < 7r/4 and are 
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denned as follows: 



NM = |i>|i> 



(60) 
(61) 
(62) 
(63) 
(64) 



|V>2> 
|V>4> 
|V>6> 
1^8) 



|0) (cos 6>i|0) +sin0i|l)), 
|2)(cos 6» 2 |1) + sin# 2 |2)), 
(cos03|l)+sin0 3 |2})|O}, 
(cos<9 4 |0) + sin 6» 4 |1»|2), 



|V3} = |O)(-sin0 1 |O} + cos0 1 |l}), 
|^ 5 > = |2)(-sin0 2 |l> + cos0 2 |2», 

|V 7 ) = (-sin0 3 |l)+cos0 3 |2))|O), 
|^ 9 ) = (- s in0 4 |O) + cos0 4 |l))|2). 



Let 83(61, 62, 03,04,) denote the rotated domino basis parametrized as above. Then the original 
domino basis is S3 := S' 3 (7r/4, 7r/4, 7r/4, 7r/4). 

Reference [BDF+99] shows that states from the domino basis S3 cannot be perfectly discrim- 
inated by asymptotic LOCC and conjectures that the same holds for the rotated domino basis 
S3(0i, 02,03, 64) f° r an y < 0\, 02, 03, 04 < 7r/4. In the next section we give an alternative proof 
that quantifies the nonlocality of the original domino states S3 and then adapt the argument to 
the rotated domino states, thus resolving the conjecture. 

5.2 Nonlocality of the domino states 

To lower bound the nonlocality constant of the domino states S3, we put an upper bound on their 
rigidity. In other words, we show that measurement operators that only slightly disturb these states 
are weak (approximately proportional to the identity operator). The key ingredient of the proof is 
Lemma 5 from Section 4. 

Lemma 6. The domino state basis S3 is A-rigid. 

Proof. The claimed result can be restated as follows (see Definition 7): 



where i,j, k,t £ {0, 1, 2} and i 7^ j or k 7^ t in the second equation. First we prove the bound for 
the diagonal elements and then we proceed to bound the off-diagonal ones. 

Bounding the diagonal elements: 

We start by bounding the differences of the diagonal elements of matrices a and b separately. Let 
us rewrite the definition of 5 from Equation (18) in the case of product states = \a.i)\(3i): 




-Tr(a®b) < 4<5Tr(a <g> b), 
9 

\aijb kt \ < 45Tr(a (8) b), 



(65) 



(66) 



iHgM mm)\ 



(67) 




(68) 



aool v 7 ( & oo + b i + 610 + &11X&00 - froi - ^10 + hi) 



boo -hi + 2ilm6i 



(69) 



\/{boo + b n ) 2 - (boi + b w ) 2 



\bpo - b u 
\boo + b u 
\boo - 611 



(70) 



(71) 



- Tr(6) 
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Note that the cancellation of |aoo| is valid since aoo 7^ by the definition of stage I. Applying a 
similar argument to the pairs of states from the other three tiles of size 2, we get that for any 
i€{0,2}, 



5 Tr(a) > \an — a» 



and 



<5Tr(6) > \b 



11 — Ojj . 



(72) 



Using these bounds and the triangle inequality, we can bound the difference between the first and 
last diagonal elements: 



|«oo — 022! < |aoo — a n| + l a n — 022! < 2(5 Tr(a) 



(73) 



and similarly 1 600 — £>22 1 < 2(5 Tr (6). 

Next, we use the bounds on the differences of the diagonal elements of a and b to bound the 
differences of the diagonal elements of a (g) b. For all k,t G {0, 1,2} we have 



\aubjj — akkht\ < \<Hibjj — akkbjj\ + \a-kkbjj — akkht\ 



I «!! — flfcfc I + I «fcfc I • I bjj — bu I 



< \bjj\ ■ 2<5Tr(a) + \a kk \ ■ 2<5Tr(6) 

< 4<5Tr(a<g>6). 



(74) 
(75) 
(76) 
(77) 



Using this inequality we can obtain the desired bound (65) for the diagonal elements: for all 
i,j G {0, 1, 2} we have 



a a b n ~ \ Tr ( a ® & ) 
y 



fe,te{o,i,2} 



< 1 £ 



\aubjj — akkht\ 



fc,te{o,i,2} 
< 4(5Tr(a<g>6). 



(78) 

(79) 
(80) 



Bounding the off-diagonal elements: 

From Lemma 5 we know that y/\T\\ ■ \T 2 \ 5 Tr(a®6) > | o ri r 2 1 ~\b Cl c2 1 ) where T\ and T2 are two distinct 
tiles and \T t \ is the area of the tile containing (r t , q). For (r%, c{) = (1,1) and any (7^, C2) 7^ (1,1) 
we get 

V25Tr(a®b) > \a nr2 \ • |6 C1C2 |. (81) 
Similarly, for any (7*1, ci) and (r2,C2) that belong to distinct tiles of size two we get 

2(5Tr(a<g)6) > \a nr2 \ ■ |6 ClC2 |. (82) 

Now it only remains to bound the following four off-diagonal elements (each of which corresponds 
to one of the four tiles of size 2): 

|aoo| • |6oi|, Kil • I&22I, I a 22|-|6i2|, |ai2|-|6 o|- (83) 
To bound |aoo| ■ |6qi|, first choose (r2,C2) = (1,0) and use Equation (81): 

\/25 Tr(a <g) b) > \a n \ • |&io| = l«n| • |&oi|- (84) 
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Now it only remains to replace on by aoo- Notice from Equation (72) that 5Tr(a) > |an — aoo| > 
I aoo I - |«n|, so 

\/26Tr(a®b) > |an| • |6 i| > (|a 00 | - 8Tr(a)) ■ \b i\ (85) 

> |«oo| • |&oi| -8Tr(a®b) (86) 

where the last inequality holds since | £»oi I < max{6oO)&n} < Tr(6) as b is positive semidefinite. 
After rearranging the previous expression we obtain 



(1 + V2)<5Tr(a®o) > | aoo | • IM- (87) 

By appropriately choosing the value of (r2, C2) and using a similar argument, we get the same upper 
bound for the remaining three off-diagonal elements listed in Equation (83). Since the constants 
obtained in bounds (81), (82), and (87) satisfy max {V2, 2, 1 + \/2} < 4, we have shown that 
Equation (66) holds for all off-diagonal elements of a b. □ 

Together with Equation (35) this implies that the nonlocality constant for the domino states is 
77 > 1/8. To get an explicit value for the lower bound on the error probability, we use Theorem 3 
with n = 9, L = 2, and c = 4. 

Corollary 1. Any LOCC measurement for discriminating the domino states S3 errs with probability 

Pevror > 1-9 X 1(T 8 . (88) 

5.3 Nonlocality of irreducible domino-type tilings 

Lemma 6 can be easily generalized to product bases that are similar to domino states on larger 
quantum systems. 

Lemma 7. LetdA,ds > 3 and let S be an orthonormal product basis of C dA ®C dB . If S induces an 
irreducible domino-type tiling of diameter D then S is 2D-rigid (see Section 4-1 for terminology). 

The proof is similar to that of Lemma 6 and appears in Appendix A. 

To bound the error probability, we use Theorem 3 with n = d^ds, L = 2, and c = 2D. 

Corollary 2. Any LOCC measurement for discriminating states from an orthonormal product basis 
of C dA C rfs that induces an irreducible domino-type tiling of diameter D errs with probability 

1 

Pcrror " 216Di(d A d B )5- {89) 

5.4 Nonlocality of the rotated domino states 

The following is an analog of Lemma 6 for rotated domino states. 
Lemma 8. The rotated domino basis S3(#i, #2, #3) #4) * s s in26> -r %^ where 



C := 6(l + 6V2 + 2y 3(6 + y/2)j < 114 (90) 

and 8 := min {6q, 62,83, 64}- 
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The proof appears in Appendix B. 

Again, we use Theorem 3 to lower bound the error probability. Here the parameters are n = 9, 
L = 2, and c= 114/sin(20). 

Corollary 3. Any LOCC measurement for discriminating S%{9i, 9i, #3, #4), the set of rotated 
domino states, errs with probability 

Perror > 2.4 x KT 11 sin 2 (26l), (91) 

where 9 := min {9\, 02, #3, 9 a}. 

Note that as 9 approaches zero, the rigidity bound tends to infinity and the bound on the error 
probability goes to zero. As the original domino basis is transformed continuously to the standard 
basis, the nonlocality decreases to zero. Moreover, since any orthonormal product basis of C 3 <S> C 3 
is equivalent to Ss(9i, 92, #3, #4) (up to local unitary transformations) for some angles 9i [FS09], 
Corollary 3 effectively covers all product bases of C 3 ® C 3 . 



6 Limitations of the framework 

6.1 Dependence of the nonlocality constant on n 



Recall that in Theorem 2 we established the lower bound p err0 v > 011 ^ ne error probability, 



_r 2 

27 1 

where 7/ is the nonlocality constant and n is the number of states. Intuitively it seems that it should 
be possible to prove a stronger lower bound on p e rror as n increases. However, to lower bound p e rror 
by a fixed constant in any dimension using our framework, one would have to prove a lower bound 
on 77 that increases with n. 

Let us consider the problem of discriminating orthonormal product states. In the next lemma 
we show that it is not possible to obtain such strong error bounds using our framework in its present 
form. We do this by proving a fixed upper bound on the nonlocality constant in any dimension. 

Lemma 9. Let S be a set of orthonormal product states in C d - A (g> C dB . The nonlocality constant 
of S satisfies 77 < 2. 

Proof. Let n = \S\ and \ipi) = |aj)|/3j). Fix some small e > 0, choose any i E [n], and define 

a = \ai)(ai\+eI dA , b = + el da . (92) 

Note that a and b have full rank and are positive semidefinite. We can easily check that 

Tr(a) = l + ed A , Tr(6) = 1 + ed B , max(^ fc |(o ® b)\^ k ) = (1 + e) 2 . (93) 

ke[n] 

Using these observations together with the definition of n in Equation (19), we get 

/ (1 + e) 2 1\ / max fc g [w] (-0 fc |(q® b)\ip k ) _ 1 \ , , 

\(l + ed A )(l + ed B ) n)-\ E ieW (^l(« ® b)\^) n) [ ] 

= 5 s (a®b) (95) 

< 1, (96) 

where the last inequality follows directly from Definition 5. As e — > 0, the left-hand side goes to 
7/(1 — ^). We can choose e arbitrarily small, so 7/(1 — ^) < 1 and thus 7/ < = 1 + < 2 since 
n > 2. ' " □ 
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6.2 Comparison to the result of Kleinmann, Kampermann, and Brufi 

The main application of the framework introduced in this paper is to show the impossibility of 
asymptotically discriminating a set of states S with LOCC. We do this by showing that the non- 
locality constant of S is strictly positive. In other words, the nonlocality constant being zero is a 
necessary condition for the sates in S to be asymptotically distinguishable with LOCC. Another 
necessary condition is presented in recent work of Kleinmann, Kampermann, and Brufi: 

Theorem ([KKB11]). Let S = {pi, . . . ,p n } be a set of states such that f^kerpj does not contain 
any nonzero product vector. Then S can be asymptotically discriminated with LOCC only if for 
all x with < x < 1 there exists a positive semidefinite product operator E satisfying all of the 
following: 

1. E i Tr(£? /0i ) = l, 

2. maxj Ti(Epi) = x, 

3. Tr(EpiEpj) = for any i ^ j. 

One should note however that in contrast to the above qualitative result, our framework can be 
applied to any set of orthogonal pure states (with no restriction on p| i ker pi) and can be used to 
obtain explicit lower bounds on the error probability. It is an open question whether our necessary 
condition ("the nonlocality constant of S is zero") or that of the above theorem is also sufficient. 
The lemma below shows that if our necessary condition is also sufficient then so is that of [KKB11]. 

Lemma 10. Let S = {iV'Oligfn] be a se ^ of orthogonal pure states such that ker (\ipi){ipi |) does not 
contain any nonzero product vector. If for all x with 1/ra < X < 1 there exists a positive semidefinite 
product operator E satisfying conditions 1-3 from the above theorem, then the nonlocality constant 
r] of S is zero. 

Proof. Consider x € (-, rrzr) an d a positive semidefinite product operator E x satisfying conditions 
1-3. Conditions 1 and 2 imply that {ipi\E x \ipi) > thus making 5s(E x ) well defined (see Defini- 
tion 5). Moreover, by condition 3 we have that \(^i\E x \ipj}\ 2 = for all i ^ j. Hence 5s(E x ) = 
according to Definition 5. Finally, from conditions 1 and 2, we get that 

m&Xi(ipi\E x \ipi) _ maxjTr(E pi) _ 

Using these observations we can rewrite Equation (19) in the definition of rj as 

n(x-^)<0. (98) 
Since x > \ h follows from the above inequality that rj = 0. □ 

7 Discussion and open problems 

We have developed a framework for quantifying the hardness of distinguishing sets of bipartite 
pure states with LOCC. Using this framework, we proved lower bounds on the error probability of 
distinguishing several sets of states, as summarized in Table 1. 

This work raises several open problems. While we were able to lower bound the nonlocality 
constant rj in many cases, it could be useful to develop more generic approaches to computing or 
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Set of states 


c 


V 


Perror 


Domino states 


4 


1 

8 


1.96 x 1(T 8 


Domino-type states 


2D 


1 

Id 


1 

2\QD 2 {d A d B f 


^-rotated domino states 


114 


sin 29 


2.41 x lO -11 sin 2 (20) 


sin 29 


227 



Table 1: Rigidity c and lower bounds on the nonlocality constant 77 and error probability p e rror f° r various 
states. 



lower bounding this quantity. We are also interested in applying our method to other sets of states. 
For example, we would like to apply the method when S is an incomplete orthonormal set (e.g., 
the domino basis with the middle tile omitted) or a product basis with tiles of size larger than 
two (see Fig. 7 for concrete examples of such tilings where no bounds on p e rror are known). It is 
unknown whether there exists a set S of 2-qubit states that can be perfectly discriminated with 
separable operations, but for which any LOCC protocol has Perror(S) > (see [DFXY09] for all 
possible candidate sets). Finally, it would be interesting to consider random product bases, since 
this would tell us how generic the phenomenon of nonlocality without entanglement is. 



□ 




□□ 
□ 

L ) 






□ 




□ 


■n 


□ 



Figure 7: Tilings corresponding to an incomplete orthonormal set in C 3 ® C 4 (left) and a product basis of 
C 5 ® C 5 with larger tiles (right). On the right, the tiles of size four are induced by states of the form |±)|±) 
and one of the tiles corresponds to the four corners of the grid. 

We discussed some limitations of our framework in Section 6, but we would like to better 
understand how broadly the framework can be applied. In particular, can it always be used to 
obtain a lower bound on terror whenever such a bound exists? For example, from Section 5.4 we 
know that the answer to this question is "yes" for orthonormal product bases on two qutrits. 

Finally, the gaps between the classes of separable and LOCC operations exhibited by our frame- 
work are rather small (see Table 1). One cannot hope to do significantly better within our frame- 
work, as shown in Section 6.1. Is this due to limitations of our framework or because orthonormal 
product states in general can be discriminated well by LOCC? 

Along these lines, a major open question raised by our work is the following: does there exist 
a sequence Si, S2, S3, . . . of sets of orthonormal product states such that 

lim rf£g c (Si) = 1? 

Existence of such a sequence would give a strong separation between the classes of separable and 
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LOCC measurements. Note that the local standard basis measurement followed by guessing gives 
the correct answer with probability at least 1/Li, where Li is the maximum number of states within 
a tile in the tiling induced by Si. Thus for any such sequence, the value of L\ must grow with I. In 
particular, the number of states (and hence the local dimensions) must also grow with /. 
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A Rigidity of domino-type states (Lemma 7) 

Lemma 7. LetdA,d,B > 3 and let S be an orthonormal product basis ofC dAl S>C dB . If S induces an 
irreducible domino-type tiling of diameter D then S is 2D-rigid (see Section J^.l for terminology). 

Proof. We mimic the proof of Lemma 6 and make the appropriate generalizations when necessary. 
We want to show that 



a a b jj ~ T~T Tr ( a ® b) 
dAUB 



<2D5Tr(a®b), (99) 
\<kjbkt\ < 2D6Tr(a(g)b), (100) 

where % ^ j or k ^ t in the second inequality. 
Bounding the diagonal elements: 

Using the calculation in Equations (68-71) we can bound the difference of diagonal entries of a and 
b. Whenever there is a 2 x 1 tile that connects rows i and j, we get that 

\aa — cijj\ < <5Tr(a). (101) 

A similar equation holds for b whenever there is a 1 x 2 tile that connects columns i and j. 

Since T is irreducible, the row graph of T is connected. Moreover, any two vertices of this graph 
are connected by a path of length at most D. We apply the triangle inequality along this path in 
the same way as in Equation (73). After at most D — 1 repetitions we get that for any i and j, 



l 13 



<D8Tr(a). (102) 



A similar equation holds for b. When we repeat the calculation in Equations (74-77), we get that 
for any i,j, k, t, 

\aubjj - akkht\ < 2D5Ti(a (g> b). (103) 

Finally, we repeat the calculation in Equations (78-80) and get the desired bound stated in Equa- 
tion (99). 

Bounding the off- diagonal elements: 
From Lemma 5 we get that 

2STr(a ®b)> \a rir2 \ ■ |6 C1C2 | (104) 
for all {r\,ci) / (?"2,C2), except when (ri,c%) and (r2,C2) belong to the same tile of size two. 
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Suppose that we want to bound \a rr \ ■ \b ClC2 \ where (r, c±) and (r, 02) belong to the same 1x2 
tile. Since T is irreducible, we can find a row r' such that (r',ci) and (r',C2) belong to different 
tiles (if {r'} x {01,02} is a tile for each r' then {01,02} is a connected component of the column 
graph of T, contradicting the irreducibility of T). From Equation (82) we get 

2<5Tr(a<g> b) > \a r > r >\ • |6 C1C2 |. (105) 

According to Equation (102), L>c)Tr(a) > \a rr — a r i r >\ > \a rr \ — |a r v|. Using this observation we 
repeat the calculation in Equations (85-86) and obtain 

26Tr(a®b) > \a rr \ ■ \b ClC2 \ - D6Tr(a®b). (106) 

After rearranging terms we get 

(D + 2)5 Tr(o (g> b) > \a rr \ ■ \b ClC2 \. (107) 

The same bound also holds for entries corresponding to 2 x 1 tiles. Together with Equation (104) 
this establishes the desired bound in Equation (100). □ 



B Rigidity of rotated domino states (Lemma 8) 

In this section we prove an analog of Lemma 6 for rotated domino states 83(61, 62 , #3, #4). For 
simplicity we consider only the set 83(6) := 83(9,6,6,6) and obtain a bound as a function of 6. In 
the more general case one can choose 9 := min {#1, 9<i, ^3, #4} and use the same bound. 



Lemma 11. For j G {0,2} we have 



til 



nl ~ sin 29 



(^ll&IL + lReM)- 



(108) 



The same inequality holds for a. 



Proof. We show how to get the bound on b for j = 0. The remaining three cases are similar. 
We use the states \ip2) and ^3) from Equation (61) in the definition of 5 in Equation (67): 



[cos 9 sin 9) 



boo b i\ ( — sin I 
&10 bu) \ cos6> 

2 a 1 l _2, 



> 



(611 — 600) sin 9 cos 9 — b\o sin 9 + 601 cos 6\ 
(611 — 600) sin 9 cos 9 + Re 601 (cos 2 — sin 2 0) + i Im 6, 
&11 - &00 



01 



sin 29 + Re 6 i cos 29 



> |6ll ~ boo| sin2fl-|Re6 i|. 
By rearranging terms we get the desired bound. 

Lemma 12. If an > - \\a\\ for some s > then for j E {0,2} we have 



\b n \ < V2s5\ 



< 2(1 + V2s) 



T/ie same statement holds when the roles of a and b are exchanged. 



(109) 

(110) 

(111) 
(112) 

(113) 

(114) 
□ 

(115) 
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Proof. We show how to get bounds on b for j = 0. The remaining three cases are identical, except 
one has to use states from different tiles. 

We use Lemma 5 with tiles corresponding to states and \ipi): 

V25 \\a <8> 611^ > |an6oi| > - NL I&01I, (116) 

where the second inequality follows from our assumption \a\\\ > - Wa]]^. By rewriting this we get 
the first bound: 

l&oil^V^H&IL- (117) 
Since |Re &01 1 < IM < V%sS ||6|| , we get the second bound from Lemma 11. □ 

Lemma 13. If an > \ \\a\\ and b\\ > \ for some s > then 



I 



Tr(a ®b) 9 



Proof. We follow the proof of Lemma 6 and show the following generalizations of Equations (65) 
and (66): 



<8(1 + V2s)— 1^(8)611^, (119) 



\a>ijbkt\ < max {v^,2 J V2*}«y||o®6|| 0O . (120) 



Note that the second inequality is stronger than we need, since 1/ sin 29 > 1. 

First, we use Lemma 12 to upper bound the difference of diagonal entries of a and b. We use 
these bounds in the same way as in Lemma 6 to upper bound the differences of diagonal entries of 
a<8> b and to get Equation (119). Finally, we use Lemma 5 to upper bound most of the off-diagonal 
entries of a ® b and Lemma 12 to upper bound the remaining ones. This gives us Equation (120). 

Bounding the diagonal elements: 

From Lemma 12 we get bounds on \b\\ — ba\ and \a\\ — aa\ for i £ {0,2}. Using the triangle 
inequality, we get 

\au-a j3 \ < 4(1 + ^)— -fat (121) 

for any i, j G {0, 1, 2} (and the same for b). Using the triangle inequality once more we can bound 
the difference of any two diagonal entries of a (8> b: 

\aubjj - a kk b tt \ < 8(1 + y/2s)— — \\a 6^ . (122) 

sin 29 

From this we obtain Equation (119) in the same way as in Lemma 6. 
Bounding the off- diagonal elements: 

Equation (120) can be obtained from Lemma 5. For most of the entries the constant is either \[2 
or 2, depending on the sizes of the tiles. For the remaining four entries, listed in Equation (83), we 
proceed in a slightly different way. For example, for aoo^oi we use Equation (117) to see that 

KoHM < \H^-^sb\b\\^. (123) 

A similar strategy works for the remaining three entries. □ 
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Lemma 14. Fix any s > 3 and let 

r(s) \U\3 S ;'2(l + ^/2 S )V3 sjf { ' 

V shfW ^ F£J then °H ^ 7 H a lloo arid &n > 7 IHIoc" 

Proof. We get one of the two lower bounds almost for free. We combine this with Lemma 5 and 
the triangle inequality to get the other lower bound. 

If maxj an = an then an > | Tr(a) > | Halloo > \ 1 1 a. 1 1 00 and we are done with a. Similarly, 
if max, bu = bn then bu > ^ H&Hoo- Thus it only remains to consider the cases when maxj an E 
{aoo, 0,22} and maxj bn S {boo, ^22}- By symmetry, it suffices to consider the case where maxj an = 
a22 and maxj bn = boo- The remaining three cases are similar. 

Using the tiles that correspond to states \1p6j) and 1^4,5)) we get 

25 \\a ® 611^ > |a 22 & i| > - Ha^ |6 i|. (125) 
Thus |Rei>oi| < l^oil < 6£ ||6|| and using Lemma 11, we get 



boo - hi < \hi - b o\ (126) 
2 

sin W" ""' x 

L-l 
sin 2(9 



< -^(^ ll&lloo + |Re6 01 |) (127) 

< 14— 4t77 l^lloo - (128) 



We assumed that maxj bn = boo-, so 



l\\b\\ 00 <boo<b 11 + U^\\b\\ 00 . (129) 

By assumption, -J^e < f^j < jg (§ ~ j) , so we get the desired bound &n > i Ho^. 
As we have a lower bound on bn, we can use Lemma 12 and get 

hi-a 22 | <2(l + v / 2s)— -wlklL- (130) 

sin 20 00 

We assumed that maxj an = a>22 , so we can rewrite this as 

~ HalL < a 22 < a n + 2(1 + V2s)-^ . (131) 

By assumption, < < 2 (i+V2s) (I ~ *)' so we get the desired bound fl n > 7 IklU- n 

Lemma 15. For any fixed s > 3 we have the following: 

• Vliro<lhthen ^ - f < 8(1 + yfa)^, 



*/ > 4t then " i M - | < r(s) 

J sm20 — r(s) lri</ M 'I — 



Tr(a®b) 9 



'sin 29' 



where r(s) is defined in Equation (124). 

Proof. The first part follows by combining Lemmas 13 and 14. To obtain the second part, notice 
that all diagonal entries of are at most 1. Since this matrix is positive semidefinite, the 

off-diagonal entries are also at most 1, so the bound follows. □ 
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Lemma 8. The rotated domino basis S^{6\, 62,9^^64) is . c 2g -rigid where 



C := 6M. + 6V2 + 2y 3(6 + y/2)\ < 114 (90) 

and := min {#1, 62, 83, #4}. 

Proof. Let us denote the largest of the two constants in Lemma 15 by 

C(s) :=max|8(l + V2s),r(s)\ (132) 

S 8 (l + V2s),U^-,2(l + V2s)^-), (133) 
[ s — 3 s — 3 J 



: max 



where we substituted r(s) from Equation (124). We want to make this constant as small as possible, 
so the best possible value is 

C = min C{s) (134) 

s>3 

3s 

= min 2(1 + s/2s) (135) 

s>3 S — 3 



6(l + 6\/2 + 2W3(6 + y/2)), (136) 



where the minimum is reached at s = 3 + a/9 + 3/^- 6.33. □ 
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